Structural phase transitions in epitaxial perovskite films 
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Three different film systems have been systematically investigated to understand the effects of 
strain and substrate constraint on the phase transitions of perovskite films. In SrTiOs films, the 
phase transition temperature Tc was determined by monitoring the superlattice peaks associated 
with rotations of TiOg octahedra. It is found that To depends on both SrTiOa film thickness 
and SrRuOs buffer layer thickness. However, lattice parameter measurements showed no sign of 
the phase transitions, indicating that the tetragonality of the SrTiOs unit cells was no longer a 
good order parameter. This signals a change in the nature of this phase transition, the internal 
degree of freedom is decoupled from the external degree of freedom. The phase transitions occur 
even without lattice relaxation through domain formation. In NdNiOa thin films, it is found that 
the in-plane lattice parameters were clamped by the substrate, while out-of-plane lattice constant 
varied to accommodate the volume change across the phase transition. This shows that substrate 
constraint is an important parameter for epitaxial film systems, and is responsible for the suppression 
of external structural change in SrTiOa and NdNiOs films. However, in SrRuOa films we observed 
domain formation at elevated temperature through x-ray reciprocal space mapping. This indicated 
that internal strain energy within films also played an important role, and may dominate in some 
film systems. The final strain states within epitaxial films were the result of competition between 
multiple mechanisms and may not be described by a single parameter. 

PACS numbers: 68.18.Jk, 68.55.Jk, 77.55.-l-f, 77.80.-e 



I. INTRODUCTION 

Perovskite thin films have received great interest in re- 
cent years due to a variety of interesting properties such 
as high-Tc superconductivity, colossal magnetoresistiv- 
ity (OMR), ferroelectricity and metal- insulator transi- 
tion. While bulk crystals of perovskites have been well 
studied, thin films are not understood in nearly the same 
details. It is well known that the properties of thin films 
can be quite different from corresponding bulk materi- 
als. The reason is often believed to be the strain and 
high concentration of defects in the thin films. This 
work illustrates that another important parameter for 
films is the geometrical constraint imposed by the sub- 
strates. Due to the excellent epitaxy, the films are tied 
to the underlying substrate, therefore the in-plane lat- 
tice parameters of films are not free to attain their bulk 
equilibrium values. The macroscopic sizes of films are al- 
ways the same as that of substrates This is known as the 
substrate clamping effect, or substrate constraint, which 
forces the temperature dependence of in-plane lattice pa- 
rameters of films to follow that of substrates. This effect 
is often neglected but as discussed in this paper, it has 
great effect on phase transitions and domain formation 



in films. The strain and substrate constraint give us a 
method to fine-tune the properties of the thin films, or 
even create a new phase that is not possible in bulk form. 
In recent years there are emerging theoretical efforts to 
treat the effects of strain^iiSi^ based on thermodynamic 
analysis using Landau theory of phase transition. 

SrTiOs (STO) has long been an important model sys- 
tem for condensed matter physics. More recently it has 
received attention for having a large and variable di- 
electric constant, making it ideal for tunable microwave 
devices.^ Bulk STO crystals are cubic perovskite at room 
temperature, with space group Pm3m (O^), but become 
tetragonal, space group I4/mcm (D^^), below 105K. This 
phase transition was extensively studied in the 1960's and 
1970's as the prototypical example of a soft-mode phase 
transitionji^ involving a progressive softening of the zone 
boundary phonon mode due to the TiOg octahedra ro- 
tation about a formerly cubic axis. When the energy 
for this mode at the zone boundary reaches zero, the 
lower symmetry structure condenses. Since the phonon 
condensation occurs at the zone boundary, this phase 
transition is non-polar and antifcrrodistortive. The high 
temperature structure of SrTiOa is the same as that of 
BaTiOa (BTO), the canonical ferroelectric crystal. In 
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BTO, the ferroelectric phase transition is associated with 
a softening of a zone-center phonon mode in which the 
central Ti ion shifts with respect to the oxygen cage. 
For SrTiOs, the ferroelectric phase, which almost oc- 
curs at lower temperatures, is suppressed by quantum 
fluctuationsii However the ferroelectric state can be in- 
duced under pressure,^ applied electric fields,^ by doping 
impurities,^'' isotope exchange^ or in thin films, in- 
dicating that this phase transition is sensitive to lattice 
strain. Therefore, studying the properties of perovskite 
films under strain conditions may as well help to better 
understand the ferroelectric phase transitions in related 
systems, such as BaTiOs and PbTiOa. In fact, simi- 
lar structural phase transitions occur in most perovskite- 
based oxide crystals, for example SrRuOa. Therefore, 
the results we obtained here may be generalized to other 
perovskite films systems. 

SrRuOs (SRO) is a ferromagnetic perovskite with 
Curie temperature of about 160K. It has excellent elec- 
trical conductivity and is chemically stable. Since SRO 
has good lattice matching with other perovskites, it is 
the material of choice for electrodes of perovskite-based 
devices. However, a bottom SRO layer may change the 
strain state and thus the properties of the top layer. A 
detailed study on effects of buffer layers is necessary, so 
that we can either minimize the negative influence, or uti- 
lize the effect to tune the strain state of the layers above. 
Here we investigated the influence of SRO buffer layers 
on the phase transition of STO films, and the domain 
structures of SRO layers. 

We also studied the behavior of NdNiOa (NNO) thin 
films under biaxial strain. Perovskite nickelates (RNiOa 
with R = rare earth) show metal-insulator (MI) tran- 
sitions with the transition temperature rising systemati- 
cally as the size of the rare earth decreases, which implies 
an increase in the distortion away from the ideal cubic 
perovskite. This transition is first-order in nature and 
associated with volume change across the phase transi- 
tion. Bulk NNO has shown large pressure dependence 
of the MI transitioufii with change in Tc of over lOOK 
under high pressure. Similar effects have been observed 
in epitaxial NNO filmsji^ but the strain state across the 
MI transition and the microstructure of NNO films are 
not clear. 

In earlier work, we reported experimental results of the 
strain effects on the phase transitions of STO thin films, 
and the observations of the smooth evolution of the out- 
of-plane lattice parameters across the transitionii^ A the- 
oretical analysis has also been undertaken to understand 
the strain relaxation mechanism in epitaxial STO films 
and the infiuence of the SRO buffer layersi^ In this pa- 
per, a systematic study of the misfit strain and substrate 
clamping effects on phase transitions of perovskite films is 
reported. We show that the thickness of both thin films 
and buffer layers may influence Tc, and the tempera- 
ture dependence of the lattice relaxation is determined 
by the substrate clamping effect and the volume change 
associated with the phase transition. The mechanisms 



involved in determining the final strain states in films 
are discussed. 



II. EXPERIMENT 

The STO films were grown by pulsed laser deposition 
(PLD) on LaAlOa (LAO) single crystal substrates. Be- 
fore growing the STO films, SrRuOa buffer layers of var- 
ious thickness were deposited. The SRO buffer serves as 
both a method to tune the strain state in the STO films 
and as a base electrode for electrical measurements. The 
samples are of the STO/SRO/LAO type, with STO film 
thickness range from 10 nm to 1 /Ltm and SRO buffer lay- 
ers from nm to 350 nm. Details of the sample prepara- 
tion are available elsewhere. We report on experiments 
of both the STO and SRO layers. 

The NdNiOa thin films, with thickness of 60 nm and 
200 nm, were grown on LAO substrates by PLD. A KrF 
excimer laser was used with an energy density of ^2.0 
J/cm^ and a repetition rate of 5 Hz. The substrate was 
heated to 700°C and an oxygen pressure of 200 mTorr was 
used during the deposition. After deposition, the films 
were cooled to room temperature at a rate of 60°C per 
minute. Resistivity was measured by the standard four 
probe method through 5 K to 300 K with current about 
0.5 ^A. 

X-ray diffraction (XRD) measurements were carried 
out at beamline X22A and X22C at the National Syn- 
chrotron Light Source, Brookhaven National Laboratory. 
Synchrotron radiation has the advantages of high inten- 
sity and small divergence angle, which are critical for 
studying very thin films grown on substrates and buffers 
of similar structures. X22A has a bent Si(lll) monochro- 
mator, giving a small beam spot and a fixed incident pho- 
ton energy of 10 keV. The longitudinal resolution with a 
Si(lll) analyzer was at least 0.001 A^^ (HWHM) for 
an (002) peak, as measured from the substrate. X22C 
has a Ge(lll) double monochromator, focused in ver- 
tical plane, with variable energy range 3^12 keV. The 
samples ware cooled in a closed cycle refrigerator with 
a temperature control better than ±0.5 K. Above room 
temperature the sample was heated in a high tempera- 
ture capable displex with a base temperature near 10 K 
and a maximum of 800 K. 

Throughout this paper we list most peaks with refer- 
ence to a cubic cell, unless otherwise specified. For cubic 
notations, we define the axis normal to the surface of the 
substrate is the c axis. Thus our scattering peaks, and 
reciprocal lattice positions, are listed as {h k I) with I 
the reciprocal direction perpendicular to the surface. For 
tetragonal or orthorhombic phases, a subscript t or o is 
added to the peak index, and the orientation is redefined 
accordingly. 
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FIG. 1: Raw XRD scan for (0 2) peaks on 200 nm STO 
/ 350 nm SRO / LAO sample, showing well resolved peaks. 
Photon energy is 10 keV. 
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FIG. 2: Temperature dependence of (1/2, 1/2, 7/2) STO su- 
perlattice peak for 1 ^m STO / 350 nm SRO / LAO sample. 
The other peak is from SRO layer. 



III. RESULTS 



SrTiOp 



3 



The high intensity and the coUimation of the syn- 
chrotron x-ray source allows for scattering from very thin 
films in a high resolution mode to easily separate film- 
buffer-substrate peaks. A typical scan over the (0 2) 
peaks is shown in Fig. ^ The STO and SRO peaks are 
well resolved, moreover the bulk LAO peak can be used 
as an internal reference to minimize systematic error. X- 
ray diffraction and TEM (not shown here) show that all 
the thin films are in excellent epitaxy with the LAO sub- 
strates and SRO buffer layers. 



1. The Phase Transition 

In bulk STO, the soft mode phase transition occurs 
at 105 K. This structural phase transition from cubic to 
tetragonal symmetry involves the rotation of the TiOg 
octahedra around [0 1] axis. Because the TiOg octahe- 
dra share the corner oxygen atoms, the neighboring TiOg 
rotate in opposite directions. This leads to a tetragonal 
lattice with space group I4/mcm. The unit cell of the 
tetragonal phase is about ^/2a x ^/2a x 2a, where a is 
the lattice parameter of original cubic unit cell, therefore 
the volume of the tetragonal unit cell is about four times 
as that of the cubic unit cell. This tetragonal phase is 
a superlattice of the original cubic phase. The selection 
rule of the diffraction peaks is (in tetragonal notation): 
hkl : h + k + l =_2n; Okl : I = 2n; hOl : I = 2n; 
hhl : I — 2n; hhl : I = 2n 

In pseudo-cubic notation, the additional superlattice 
peaks are at half integer positions, such as (1/2, 1/2, 
3/2),. 

In STO films, we also observed the appearance of the 
superlattice peaks below Tp, as shown in Fig. [3 which 
indicates that the phase transition also occurs in thin 
films. The temperature dependence of the peak inten- 
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FIG. 3: Temperature dependence of superlattice peak inten- 
sity on STO film thickness. Bulk value from Ref. ifoL 



sity, as seen in Fig.|31 shows obvious rounded tails around 
Tc , which means the phase transitions occur over a range 
of temperature. This is considered normal for films due 
to strain and defects, but it makes determination of Tc 
more difficult. We chose to linearly extrapolate the part 
of curves between 10% and 40% of maximum intensity, 
and define the Tc as where this line crosses the 0% line. 
This method is justified by both our data as well as in 
literaturepiS Though it might give a Tc slightly higher 
than the true transition, we believe the comparison be- 
tween different films is accurate. 
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Our first set of samples have STO film thickness vary- 
ing from 50 nm to 1 /xm, with the same SRO buffer layer 
thickness of 350 nm. This shows us how the phase tran- 
sition depends on the thickness of STO films, as other 
conditions keep the same. The Tc difference, ATc, is 
about 14 K between the 50 nm sample and 1000 nm sam- 
ple, as seen in Fig. 13 The thinner film has higher Tc 
while thicker film has Tc closer to the bulk value. The 
in-plane tensile strain measured by XRD at room tem- 
perature increases from 0.01% in 1000 nm film to 0.23% 
in 50 nm film. Therefore in this case, Tq increases mono- 
tonically with increasing in-plane strain. 

The second set of samples have the same STO thickness 
of 200 nm, but SRO thickness ranges from nm to 350 
nm. This gives us information about the influence of the 
SRO buffer layers. We found that varying the SRO thick- 
ness has even bigger effect on Tc, as shown in Fig^ The 
thinnest (5 nm) SRO layer causes the highest Tc in STO 
films, which is about 50 K above the bulk value. When 
SRO thickness increases, Tc decreases. However, room 
temperature XRD measurement shows that the in-plane 
strains in these samples are still quite small, with the 
largest one only about 0.25%, so this amount of change 
in Tc is fairly remarkable. It is interesting to notice that 
with SRO thickness decreased to nm, i.e., STO films 
grown directly on LAO substrate, the Tc goes back to 
about 130 K. This signals that the SRO layers have a 
significant effect on the strain states of STO films. 
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FIG. 4: Temperature dependence of superlattice peak inten- 
sity on samples with different SRO thickness. Bulk value from 
Ref.[l3 



2. Lattice Parameters 

The lattice constants were measured by conventional 
diffraction as well as grazing incidence X-ray diffraction 
(GIXD) technique. Due to the high intensity of the syn- 
chrotron radiation, peaks from all three layers, even as 
thin as 20 nm, can be seen clearly. Measuring multi- 
ple peaks in (0 /) direction gives a very accurate mea- 
surement of the c lattice constant. In-plane lattice con- 
stants are obtained either by measuring a partially in- 
plane peak such as (2 2)/(0 2 2) and then triangulating 
using the (0 I) values, or by measuring a (2 0)/(0 2 
0) type reflection using grazing incidence diffraction. 

We have measured the temperature dependence of 
both out-of-plane c of all samples and in-plane lattice pa- 
rameters of some samples. Fig. |5l shows the comparison 
between the lattice constants of films and STO bulk, as 
well as the LAO substrate. We believe the difference in (2 
0) and (0 2 0) is due to experimental error. It is imme- 
diately clear that the lattice parameters evolve smoothly 
from low temperature up to room temperature, without 
showing any indication of the structural phase transi- 
tion. The in-plane lattice constants follow that of the 
LAO substrate, which implies that the substrate plays 
the controlling role here. 

We attribute the smooth variation of the lattice pa- 
rameters to the constraint applied by the epitaxy of the 
film on the far thicker substrate. This substrate clamp- 



ing effect forces the in-plane lattice parameters of films 
to follow that of substrates. Therefore, even above the 
phase transition temperature, as defined by intensity of 
superlattice peaks, STO still has different in-plane and 
out-of plane lattice constants. In fact, the symmetry of 
the films at high temperature is no longer cubic, but high- 
symmetry tetragonal, possibly P4/mmm. In other words, 
now the phase transition is no longer from cubic Pm3m 
to tetragonal I4/mcm, but from high symmetry tetrag- 
onal to lower symmetry tetragonal. In bulk STO, the 
splitting of the lattice parameters at Tc results in do- 
main formation and thus lattice relaxation. However, in 
films, there is no change in the shape of the unit cell dur- 
ing the transition. Therefore, externally, the films always 
appear as in single domain, and the phase transitions oc- 
cur without lattice relaxation. Obviously the distorted 
unit cell favors the TiOg rotation, so the starting point 
of the rotation, Tc, is much higher in films. Origins of 
this effect are discussed more fully in Section Hvl 



B. NdNiOg 

Bulk NdNiOa is different from bulk SrTiOs in that it 
shows discontinuity in temperature dependence of unit 
cell volume and lattice parameters across the first-order 
MI transitioniSi To understand how they behave in films 
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FIG. 5: Temperature dependence of superlattice peak in- 
tensity, in-plane and out-of-plane lattice parameters of STO 
films. Dashed lines indicate Tc- LAO and STO bulk datai^a 
are shown for comparison. 



under the influence of strong substrate clamping effect, 
we have measured the temperature dependence of in- 
plane and out-of-plane lattice parameters, as well as re- 
sistivity, of NNO films. 

As shown in Fig. |^ for the 60 nm sample, the in-plane 
lattice parameters show the same behavior as in STO 
films, following closely with the LAO substrates. But 
the out-of-plane lattice parameters exhibits a jump at 
the MI transition, although LAO shows nothing anoma- 
lous around this range. From Fig.jTJ we can see that the 
discontinuous volume change occurs at about 140 K in 
60 nm films during warming up, in contrast to 200 K in 
bulk. This is in good agreement with the Tc obtained 
by resistivity measurements (Fig. (S)). The 200 nm sam- 
ple behaves similarly, and its Tc is about 135 K during 
cooling down. Another important feature in films is that 
the volume expansion with respect to the volume on the 
metallic side is only 0.053% for 60 nm sample, consider- 
ably smaller than the bulk value of 0.23%. The thermal 
expansion coefficient a behaves similar in bulk and films. 

In bulk NNO, the changes in all three lattice constants 
contribute to the overall volume discontinuity, while in 
epitaxial films, this is no longer feasible, since the in- 
plane lattice parameters are constrained by the underly- 
ing substrates. Assuming a volume change is necessary 
for MI transition, then in thin films, the only possible 
way is to let out-of-plane lattice constant c change ac- 
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FIG. 6: Out-of-plane and in-plane lattice parameters of 
NdNiOs films and LAO substrates during heating up, shown 
in pseudo-cubic notation. 
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cordingly. Therefore, in phase transitions with volume 
discontinuity, the substrate clamping effect still controls 
the behavior of in-plane lattice parameters, while the out- 
of-plane parameters adapt to accommodate the volume 
change necessary for the phase transition. 



C. SrRuOg 

Bulk SrRuOs undergoes phase transitions at several 
temperatures. Above 950 K (677°C), SRO is cubic with 
Pm3m symmetry. Between 950 K to 820 K it is tetrag- 
onal I4/mcm. Below 820 K, the unit cell becomes or- 
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thorhombic with space group Pbnm.22 At about 160 K, 
SRO goes through another transition to a ferromagnetic 
phase. 

We found domain structure in the SrRuOa buffer lay- 
ers at room temperature, presumably due to tetragonal- 
orthorhombic transition. This appears to be a counter- 
case to the clamping effect we mentioned before. Notice 
that we will use orthorhombic notation for SRO, but still 
employ the pseudo-cubic notation for LAO. SRO has a 
orthorhombic unit cell of about \/2axy^ax2a, where a 
is pseudo-cubic axis. 

The growth temperature 720°C (^993 K) is substan- 
tially higher than the first phase transition temperature 
of SRO. Therefore at growth temperature, SRO is pos- 
sibly cubic, with lattice mismatch fully relaxed by form- 
ing misfit dislocation. Upon cooling down to room tem- 
perature, it may form poly-domain structure to relax 
the straini22i24 According to the lattice parameters of 
pseudo-cubic LAO and orthorhombic SRO, there are six 
possible domain orientations. With long axis (c) of SRO 
unit cell parallel to the [1 0]c of LAO lattice, we have 
a pair of 90 degree twin domains, A and A'. Similarly, 
with long axis parallel to [0 1 0]c of LAO defines B and 
B' domain pair, and parallel to [0 l]c of LAO defines 
C and C" domains. 

We examined SRO layers in several samples by x- 
ray reciprocal space mapping at room temperature. At 
pseudo-cubic (113)c position, which would be (206)o for 
c domain or (422)o for a(b) domains, we observed two 
SRO peaks splitting along the cubic [110]c direction. This 
arrangement of the two peaks suggests c domains, with 
SRO c axis perpendicular to the surface of the films. The 
two peaks correspond to (206)o of one domain and (026)o 
of the other domain from a 90° twin structure. 

We have performed temperature scans from 300 K to 
700 K on a 5 nm STO / 350 nm SRO film on LAO sam- 
ple. Fig. ini Both hk plane and kl plane were scanned. 
At room temperature, there are clearly two SRO peaks. 
From 500 K the peaks began to get closer, and finally 
merged at about 600 K. Upon cooling down, the peak 
splits at the same temperature, showing a reproducible 
phase transition. This Tc is almost 150'--^200 K below the 
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FIG. 9: Reciprocal space mapping shows SRO forming do- 
main structure. The lower peaks are SRO (2 6)o and 
(0 2 6)o. The upper weak peaks are from STO (1 1 3)c. 
The unit is in reciprocal lattice unit of SRO. 
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bulk value. The temperature dependence of lattice pa- 
rameters obtained from above measurements are shown 
in Fig. EH Overall it is similar to the bulk, though the 
splitting is bigger. 
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IV. DISCUSSION 
A. Order Parameters 

According to Landau theory, the phase transitions 
of the second order can be quantitatively described by 
properly defined order parameters, which represent the 
changes in the symmetry in the crystals when they pass 
through the phase transition points. The rotational angle 
of the TiOg octahedra, has been identified as the or- 
der parameter of phase transition in ST0f2i and denotes 
the deviation from the perfect cubic perovskite structure. 
The intensity of the superlattice peak is proportional to 
the square of order parameter. For bulk STO samples, 
above the phase transition temperature Tc, the structure 
is in cubic symmetry, so there is no superlattice peak. Be- 
low T^, the TiOg octahedra rotate by an angle and 
superlattice peaks appear. Concurrently, the three lattice 
parameters are no longer equal, with c axis longer than 
the other two. The changes of the lattice parameters in- 
dicates the tetragonality of the unit cell changes at Tc- 
Therefore, we can use the tetragonality as a secondary or- 
der parameter to describe this phase transitioni ^^i^^ We 
define the tetragonality as: 

^ _ (ai-l-a2) 

^ (ai+aa+c) " 
3 

where oi, 02 and c are the lattice constants in pseudo- 
cubic notation, and the average of two in-plane values 
oi and 02 is used to account for possible experimental 
error. For the cubic phase, 7 is zero. The tetragonal 
phase has either positive or negative 7 depending on if c is 
greater than or less than a. It is interesting to notice that 
the order parameter defined by rotation angle represents 
a internal change within the unit cell, while the order 
parameter obtained from tetragonality gives the external 
change in the shape of the unit cell. 

From above results, we see an important change in 
the nature of this phase transition in thin films vs bulk. 
In bulk STO, the internal and external signals are inti- 
mately related, the TiOg octahedra start to rotate while 
the shape of the unit cell distort at Tc. Both signa- 
tures show sharp transition points at the same transition 
temperature of 105 K. In epitaxial thin films, the exter- 
nal lattice shape is no longer free to change, however, 
the TiOg octahedra still begin to rotate at particular 
temperature within this fixed cage. Therefore, although 
both the order parameters as defined before are good for 
bulk, in thin films only the internal order parameter can 
clearly indicate the phase transition. The tetragonality 
order parameter shows no phase transition, as illustrated 
in Fig. 111! In other words, the internal degree of freedom 
is now decoupled from the external degree of freedom. 

This phenomenon is not unique for the STO system. 
Similar cases happen in ferroelectric phase transitions in 
epitaxial BaTiOa or other thin films. When bulk BaTiOa 
goes through a ferroelectric phase transition at 120°C, its 



lattice constants change abruptly, while the unit cell vol- 
ume is continuous acrossing the Tq-'^'^ But in thin films, 
there is no obvious change in all three lattice param- 
eters at the phase transitioni^i*^2i2^ However from the 
dielectric measurement, Yoneda and coworkers observed 
ferroelectricity in the low-temperature phase and para- 
electricity in the high-temperature phase in films. This 
indicates that there is a phase transition with respect 
to the internal displacement of the Ti and O atoms, al- 
though the external structural component of the phase 
transition is suppressed in the presence of epitaxial strain 
and substrate constraint. 

For BTO, the primary order parameter is the sponta- 
neous polarization P (D=0), which is proportional to the 
internal relative displacement of Ti and O atoms (the fer- 
roelectricity of the material) , and again the tetragonality 
of the unit cell is the secondary order parameter 1^ We 
will see exactly the same relationship between the bulk 
and thin films as the case of STO. Therefore, from the 
view point of order parameters, we have a new kind of 
phase transition which is not possible in bulk. 



B. Domain Formation 

It is natural to ask why some films show bulk-like relax- 
ation of lattice parameters at the transitions while others 
do not. This is equivalent to asking why some films form 
domains while others do not. We give a qualitative ex- 
planation from an energetic point of view. 

If there is a material with a cubic to tetragonal phase 
transition, above Tc, cubic phase is energetically fa- 
vored, and below Tc tetragonal symmetry is favored. 
During the transition, some atoms within the material 
will have slight displacement, from XfUm positions in 
cubic phase to Xj^^^ positions in tetragonal phase. In 
other words, the potential energy minimums for these 
atoms shift, and this is the internal driving force for the 
phase transition and domain formation. If external forces 
are applied to keep these atoms stay at XfUm positions 
even below Tc, the material is strained and the total free 
energy increased. 

On the other hand, imagine we have a completely fiex- 
ible epitaxial film growth on a substrate, the film atoms 
would always align with the underlying substrate atoms. 
The substrate creates a potential energy minimum at 
each of the X^ub positions directly above the substrate 
atoms, and these minima only move with substrate. If 
we need to move some atoms away from these Xsub posi- 
tions, we must supply some energy. And this again raises 
the total free energy of the system. 

Combining these two pictures, we can see there are two 
competing energies during this process, involving internal 
strain energy from film and interface strain energy due to 
substrate. For a thick film, at the growth temperature, 
the two minimum XfUm and Xsub are aligned with each 
other through formation of misfit dislocations. When the 
temperature is below Tc, the film has the tendency to 
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FIG. 11: Order parameters for phase transitions in STO bulk and films, (a) superlattice peak intensity in bulk, (b) tetragonality 
in bulk, (c) superlattice peak intensity in films, (b) tetragonality in films. Dashed lines indicate Tc- Bulk data from Ref. Il9ll2a . 



complete the structural phase transition, and this makes 
the two potential energy minima separated into X'j:^^,^^ 
and Xsub- The deviation of an atom from one minimum 
will raise the corresponding strain energy, and we expect 
the increase to be small if the displacement is small. 

Generally the structure change is minor when temper- 
ature is near Tc, so just below Tc the potential mini- 
mum at X'j:^i^ due to film is likely to be close to Xsub 
position and possibly shallow than X^ub- Therefore the 
driving force for the domain formation is small and the 
film atoms will still follow substrate. No domain forms. 
When temperature is considerably lower than Tc and 
the depth of minimum at X'^^i^ exceeds that at Xsub, 
the stable position moves to X'^^^^, therefore forming 
domain is energetically favored. However there will be 
a barrier between the two minima, and in order to over- 
come the barrier to reach X'^^^^, the atoms need some 
thermal energy. 

For STO, there can be two scenarios, (i) since the struc- 
tural change is subtle, the potential minimum at X'^^^^ 
may always be shallower than that at Xsub] (ii) since the 
Tc is far below the growth temperature, even if the film 
has a tendency to form domain below Tc, the thermal 
energy of the film atoms is too low to overcome the bar- 
rier to the new position. So even though the minimum at 
-^'fiim is lower, the atoms still stay at the minimum Xsub 
defined by the substrate. In both cases, it appears macro- 
scopically that the film is frozen to the substrate, without 
any obvious in-plane lattice change. The domain forma- 
tion is suppressed by the substrate constraint. BaTiOa 
and NdNiOa films fall into exactly the same situations. 

For SRO, the phase transition occurs at above 820 K in 
bulk. If we assume the film has the tendency of a phase 



transition at a similar temperature, then the atoms still 
have substantial thermal energy, or mobility, to move 
around. If at a certain temperature the move from Xsub 
to Xy^i^ is energetically favorable, the film atoms will 
move to the minimum at X'^^^^^ defined by the internal 
strain energy, then new domains form. Further experi- 
ments are on going to confirm this explanation. 



C. Lattice Parameters vs Thickness 

From XRD measurements of lattice constants, we can 
see that both film thickness and buffer layer thickness 
influence the strain state of the STO films. Fig. IT^ 
and Fig. ll3lshow both the out-of-plane and in-plane lat- 
tice parameters at room temperature for different film / 
buffer combinations. For films grown on the same sub- 
strates and buffer layers, when STO film is very thin, the 
misfit strain is not fully relaxed, so the in-plane lattice 
constants are expanded, while the out-of-plane lattice 
constant is reduced in order to keep roughly the same 
unit cell volume. As the STO thickness increases, the 
strain within the STO films becomes more and more re- 
laxed, so both the out-of-plane and in-plane lattice pa- 
rameters move towards the equilibrium bulk values. The 
thickest film, the 1000 nm STO / 350 nm SRO sample, 
is already indistinguishable from cubic. Considering the 
experimental error, its three lattice parameters are al- 
most the same. Also it behaves more like bulk crystals, 
with Tc of 122K, most close to 105K. 

For the same STO films grown on different SRO lay- 
ers, the variation of STO lattice constants becomes more 
complicated. Since the lattice constants of SRO are 



9 



o 
I— 

CO 



0.393 

0.392 

0.391 

0.390 

0.389 
1.006 
1.004 
1.002 
1.000 
0.998 
0.996 
0.994 



"T- 



-T- 



o Exp. c 
□ Exp. a, 
^ Exp. 



(a) 




(b) 



Calc. In-plane a/og 

— Calc. Out-of-plane c/oj 



200 400 600 800 
STO thickness (nm) 



1000 



FIG. 12: Lattice parameters for STO films of different thick- 
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larger than that of STO, so it seems reasonable to expect 
thick SRO layers give rise to a larger strain in STO films. 
But the results turn out to be just opposite, the thinnest 
SRO layer creates biggest deformation in the STO films. 
Increasing the SRO thickness, the in-plane lattice pa- 
rameters in STO films keep decreasing. The out-of-plane 
lattice constant c is also decreasing, and always smaller 
than in-plane ones. An interesting phenomenon is all 
three parameters are larger than the STO bulk values, 
which means the volume of the unit cell is bigger than 
the bulk. 

As seen from Fig ll3l the STO lattice constants change 
drastically when SRO thickness goes from nm to about 
20 nm. If STO is grown directly on LAO substrate, the 
two in-plane lattice parameter ai and 02 are smaller than 
out-of-plane c, however, they are still larger than the bulk 
values. That is to say, STO is expanded by LAO, which 
is contradictory to the common thought. This indicates 
that the strain states in the films also rely on other fac- 
tors, such as the differences between the thermal expan- 
sion coefficients of films and substrates, and defects in 
the films. 

To understand the experiment results, we have calcu- 
lated the dependency of the STO lattice parameters upon 
the thickness of both STO and SRO layers. The details 
were published elsewhere»i& The theoretical model em- 
ployed takes into account the stress relaxation due to 
formation of orthorhombic polydomain structure in the 
SRO buffer layer, as well as the formation of misfit dislo- 
cations at the LAO / SRO and STO / SRO interfaces. It 
has been shown that the internal stress level in films can 
be controlled using buffer layer that exhibits a structural 
phase transition. 

Taking into account the multiple strain relaxation 
mechanisms, the in-plane and out-of-plane lattice param- 
eters of STO films as a function of both STO and SRO 
layer thickness can be calculated, as shown in Fig.ll2land 
Fig. 1131 In general the calculations is in good agreement 
with experiment data. 



D. Phase Diagram 

It is well known that Tc in films change with differ- 
ent film thickness or different substrates. Generally, one 
would like to parameterize the control by a simple coor- 
dinate and construct an appropriate phase diagram. The 
most commonly discussed parameter for films is strain. 
Another candidate for the STO phase transition might 
be tetragonality at room temperature. Thus we would 
like to verify a Temperature-Strain phase diagram and a 
Temperature- Tetragonality phase diagram. 

The phase transitions in strained SrTiOa epitaxial 
films have been theoretically calculated and discussed 
by Pertsev et ali The misfit strain-temperature phase 
diagram they predicted shows that, except for the high- 
temperature tetragonal phase (HT) being the distorted 
prototypic cubic phase, there are two purely structural 
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tetragonal and orthorhombic states (ST and SO) at lower 
temperature, and the phase transition temperature is 
strongly related to the misfit strain. The minimum of 
Tc is at 105K with zero misfit strain, while both posi- 
tive strain and negative strain increase Tc- 

Since there is no indication of the structural phase 
transitions with respect to the changes in lattice con- 
stants, we obtained the phase transition temperatures 
by monitoring the superlattice peaks, which corresponds 
to the rotation of the TiOe octahedra. In Fier. I14f a). we 
plot the relationship between Tp of STO films vs in-plane 
strain. The general rule for effect of in-plane strain on 
Tc is that, the larger the in-plane strain, the higher Tc. 
The strain in our samples ranges from 0.01% to about 
0.25%. In this range, the trend of our measurements are 
in general agreement with the calculations. The phase 
transition temperatures increase with increasing in-plane 
strain. Above the transition temperature, the STO film 
is a high-symmetry tetragonal phase due to the biaxial 
strain and substrate constraint. However, there are some 
differences. They predicted that Tc would increase very 
slowly for tensile strains, not higher than 115 K for 2% 
strain. But our results show remarkably larger increases 
in Tc even for much smaller strains, about 160 K for 
0.25% strain. In addition, it is possible to induce the 
same amount of in-plane strain in STO films by varying 
either STO thickness or SRO layer thickness. However 
the transition temperatures are different for these two 
cases, thus Tc vs strain is not a simple monotonic func- 
tion. 

From Fig. [TlT bl Tc vs tetragonality of STO unit cell, 
we also found there is no simple relationship between 
them. The main factor influencing the tetragonality is 
the thickness of STO itself. The thin STO films tend 
to have high tetragonality while thick STO films is more 
bulk-like with tetragonality close to 0. With the same 
STO thickness of 200 nm, the films on different SRO 
layers have almost the same tetragonality. However the 
Tcs of these films are quite different. That is to say, 
tetragonality is not an adequate factor to dictate Tc. 
This is also consistent with fact that the tetragonality of 
each film does not change across the phase transition. 

Some reasons for the discrepancy between theory and 
experimental data may include: (i) an ideal cubic sub- 
strate was used for calculations, but the real buffer- 
substrate systems usually have mosaic structures and do- 
mains. This leads to inhomogeneous strain conditions 
across the films; (ii) the different definition of Tc we 
used may give a slightly higher Tc. However, we believe 
the primary reason is because film systems are inherently 
more complicated than the model proposed. The strain 
state in a film is due to more than the difference between 
the lattice parameters of film and substrate. On one 
hand, two strain relaxation mechanisms, forming misfit 
dislocations at the film-substrate interface and forming 
domain structures within the films, greatly reduce the 
actual misfit strain comparing with that calculated by 
respective bulk values. On the other hand, the differ- 
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ence between the thermal expansion coefficients of films 
and substrates, together with substrate constraint, create 
another kind of strain, which depends on temperature, 
but not on the lattice constants of films and substrates. 
Furthermore, non-equilibrium defects increase unit cell 
volume)Si2& and thus change strain in a manner differ- 
ent than epitaxial mismatch. 

Adding together all these complications, the Tc vs 
strain is not a simple function in real films, and an equi- 
librium phase diagram may not be generally applicable. 



V. CONCLUSION 

Through a systematic study of three distinctive film 
systems, we have a more comprehensive understanding 
on the phase transitions in epitaxial perovskite films. The 
final strain state of film, and thus the Tc, results from 
competition between multiple mechanisms, and may not 
be ascribed to a single parameter, such as misfit strain. 
There are a few general rules, such as thickness depen- 
dence of strain and substrate constraint, which are fol- 
lowed by all film systems. However, a particular film 
system may possess its own unique features, which must 
be characterized individually. 

We have studied the structural phase transitions in 
strained STO films grown on LAO substrates with SRO 
buffer layers. As indicated by the superlattice peaks asso- 
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ciated with the lower temperature tetragonal phase, the 
phase transitions occur at higher temperature in strained 
films, with larger in-plane strain giving higher T^. How- 
ever, there is no obvious indication of this transition from 
temperature dependence of the lattice parameters, and 
tetragonality is no longer a good secondary order param- 
eter. This is an important change in the nature of this 
phase transitions. Due to the epitaxial strain and sub- 
strate clamping effect, the internal degree of freedom, the 
rotation of the TiOe octahedra, decouples from the ex- 
ternal degree of freedom. Our data show that both the 
thickness of STO thin films and SRO buffer layers influ- 
ence the magnitude of the strain within STO films. The 
in-plane strain decreases with increasing STO thickness, 
however, thinner SRO layers give rise to larger in-plane 
strains in STO films. This result is in good agreement 
with our calculations. Combining the Tc and strain mea- 
surements, it is found that the in-plane strains have con- 
siderably larger effect on Tc than previously predicted. 

From our results of NdNiOs films grown on LAO, it 
is much clearer that substrate constraint is an important 
parameter for epitaxial films. The in-plane lattice param- 
eters of film are tied down by the underlying substrate, 
following the trend of substrate through the whole tem- 
perature range, without change even across the phase 
transition in films. This leaves the out-of-plane lattice 
parameter the only variable to accommodate the vol- 



ume change during the phase transition. This substrate 
clamping effect has a profound influence on the strain 
within film, and can induce a substantial change in phase 
transition temperature Tc- 

SrRuOa films show us even more complications. In this 
system we observed domain formation at elevated tem- 
perature. This means the substrate constraint sometimes 
may give way to other much stronger strain relaxation 
mechanism. This phenomenon illustrates that the final 
strain states in films are the results of interaction be- 
tween multiple competing mechanisms, such as internal 
strain energy within films and interface strain energy. 
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